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Abstract 

The parton momentum density distributions of the proton are determined from 
a next-to-leading order QCD analysis of structure functions measured at HERA and 
by fixed target experiments. Also included are data on the difference of the up and 
down anti-quark distributions. The uncertainties in the parton densities, structure 
functions and related cross sections are estimated from the experimental errors and 
those on the input parameters of the fit. Several QCD predictions based on the 
parton densities obtained from this analysis are calculated and compared to data. 



1 Introduction 



Parton momentum density distributions are important ingredients in the calculation of 
high energy hadron-hadron and lepton-hadron scattering cross sections. In these calcu- 
lations the cross sections are written as a convolution of the parton densities and the 
elementary cross sections for parton-parton or lepton-parton scattering. Whereas the lat- 
ter can be perturbatively calculated in the framework of the Standard Model, the parton 
densities are non-perturbative and cannot, at present, be calculated from theory; they are 
obtained from fits to a body of high energy scattering data. In these fits measurements 
of deep inelastic lepton-hadron scattering play an important role because in this process 
the hadron structure is directly probed by the structure-less incident lepton. 

Standard sets of parton density distributions have been published by, among others, the 
groups CTEQ, MRS and GRV, see |2|, |3|] for their most recent results. These parton dis- 
tribution sets are widely used as an input to cross section calculations. However, estimates 
of the uncertainties on the parton densities are not provided. These uncertainties tend 
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to dominate the errors on the predicted cross sections and are an increasingly important 
limitation on the interpretation of hadron-hadron or lepton-hadron collider data. 

For instance with the integrated luminosity of about 50 pb _1 collected at HERA during 
the years 1994-1997 a new kinematic domain of large x and Q 2 becomes accessible for 
the study of deep inelastic scattering in ep collisions. Here Q 2 = —q 2 with q the four- 
momentum transfer from the incident electron to the proton and x = Q 2 /2p ■ q with p 
the four-momentum of the proton. Measurements by ZEUS of the e + p — > e + X neutral 
current (NC) and e + p — ► vX charged current (CC) Born cross sections have recently 
been published ||, ||] (see || for recent NC and CC cross section measurements by HI). 
These data extend the largest measured value of Q 2 from about 5000 GeV 2 [|7| to Q 2 = 5- 
6 x 10 4 GeV 2 . Standard Model predictions for the cross sections calculated with the 
parton distribution set CTEQ4 || are in good agreement with the NC results but fall 
below the ZEUS CC measurements at large x and Q 2 . This could be an indication of new 
physics beyond the Standard Model but might also be due to an imperfect knowledge 
of the parton densities in this kinematic region. For instance in || it is shown that a 
modification of the CTEQ4 down quark density yields predictions in agreement with the 
CC e + p data. 

To investigate these issues we have performed a next-to-leading order (NLO) QCD analysis 
of deep inelastic structure function data from HERA and fixed target experiments. The 
HERA data extend down to low values of x ~ 10~ 5 where charm mass effects must be 
taken into account. For this we could have used, since it is implemented in the QCD 
evolution code, the heavy quark three flavor scheme of |TIJ. However, the CC cross 
sections and the Z-exchange contribution to the NC cross sections cannot be calculated 
in this scheme. Because we wish to compare the QCD predictions to the ZEUS NC and 
CC data at large Q 2 we do not use the heavy quark scheme and assume, instead, that 
the charmed and bottom quarks are massless. As a consequence, we have restricted the 
kinematic domain of the analysis to larger values of x > 10~ 3 where the massless approach 
is expected to yield reasonable QCD extrapolations to large Q 2 , see Section |6|. 

In this paper we present the results of the QCD analysis including estimates of the un- 
certainties in the parton densities, the structure functions and the Standard Model pre- 
dictions for the NC and CC cross sections. The results of this analysis, including the full 
error information, are available from the web site given in Section |7|. 



2 QCD fit 



The data used in the fit were the proton structure function (Ff ) measured at HERA by 
| TT| together with the proton and deuteron (F^) structure functions 

Also included in the fit were measurements 



ZEUS 



and HI 



from fixed target experiments |L2], [T3|, [LJ, [15 



of the ratio F$/F% by NMC |]| and neutrino data on xFf e from CCFR JXTJ. Recent 
results from E866 |18| on di-muon (Drell-Yan) production cross sections in pp and pd 
scattering yields information on the difference between the down and up anti-quark den- 
sities. Because a NLO calculation of the Drell-Yan cross sections is not incorporated in 
the present analysis we have, instead, constrained this difference by using their results on 
x(d — u) as an input to the fit. 
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Corrections for nuclear effects were applied to _F 2 d , F^/F^ and xF^ Fe . For this purpose 
the structure function (per nucleon) of nucleus A was parameterized as 

^ = ^[1 + ^(^-1)] (1) 

where F t N is the structure function of a free nucleon (taken to be the average of proton 
and neutron) as predicted by the QCD fit, Ra is a parameterization of the ratio Ff- / 
taken from the literature (see below) and is a parameter which controls the size of 
the nuclear correction applied. 

In Fig. [L]a we plot the ratio Flf/Fg versus x as parameterized by Gomez et al |I9|| . We 



allow for an error of 100% on this nuclear correction as indicated by the shaded band 
in the figure, that is, we set K d = 1 ± 1. The nuclear correction applied to xF£ Fe is 
shown in Fig. [Tft> . For this we took the parameterization of Eskola et al [^0| of the nuclear 



correction for iron on the valence quark density x(q — q) oc xF^ . Here we allowed for an 
error of 50% i.e. Kp e = 1 ± 0.5 (shaded band in Fig. [l|b). 

The QCD predictions for the structure functions were obtained by solving the DGLAP 
evolution equations [ETJ in NLO in the MS scheme ||22|| . These equations yield the quark 
and gluon momentum densities (and thus the structure functions) at all values of Q 2 
provided they are given as functions of x at some input scale Qq. Also required is an input 
value for the strong coupling constant which, in this analysis, was set to a s {M 2 z ) = 0.118 
(±0.005). With a charm (bottom) threshold of m c n^ = 1.5 (5) GeV in the Q 2 evolution of 
a s this corresponds to values of the QCD scale parameter A(MS ) = (404, 343, 243) MeV 
for / = (3, 4, 5) flavors. 

As already mentioned in the introduction, we use in this analysis the light quark variable 
flavor number scheme where all quarks are assumed to be massless and where charm 
and bottom are dynamically generated above some given thresholds Q 2 and Q 2 . These 
thresholds were set to 4 (±1) and 30 GeV 2 respectively. 

The input scale of the DGLAP evolutions was chosen to be Qq = 4 GeV 2 and the gluon 
distribution (xg), the sea quark distribution (xS), the difference of down and up anti- 
quarks (xA) and the valence distributions (xu v , xd v ) were parameterized as 

xg(x,Q 2 ) = A g x^{l-x)^{l + lg x) 

xS{x,Q 2 ) = 2x(u + d + s) = A s x Ss (l-x) Vs (l + -f s x) 

xA(x,Ql) = x(d-u) = A a x Sa (1 -x) Va (2) 

xu v (x, Ql) = x(u — u) = A u x Su (l — x) Vu (l + ^ux) 

xd v (x,Ql) = x(d - d) = A d x 5d (l - x) m (l + j d x). 

In the fit, we assume that the strange quark distribution x(s + s) = 2xs is a given 
fraction K s = 0.20 (±0.03) of the sea at the scale Q 2 = 4 GeV 2 , in accordance with the 
measurements of CCFR |2l 



The parameters A g , A u and A d were fixed by the momentum sum rule and the valence 
quark counting rules: 

(xg + xS) dx — 1, / u v dx = 2, / d v dx = 1 
'o Jo Jo 

where xS = xS + xu v + xd v denotes the singlet quark density (sum of all quarks and 
anti-quarks). There are thus 16 free parameters describing the input parton densities. 



3 



From the evolved parton densities the relevant structure functions were calculated in NLO. 
In the calculations both the renormalization and the factorization scale were set equal to 
Q 2 . Higher twist contributions to F% and F$ were taken into account phenomenologically 
by describing these structure functions as 

F kt = f QC» [1+h{x)/q2] (3) 

where i^ CD obeys the QCD evolution equations. The function H(x) was parameterized 
as a fourth degree polynomial in x with five free parameters: 

H(x) = h + h x x + h 2 x 2 + h 3 x 3 + h 4 x 4 . (4) 



We assume in this analysis that H(x) is the same for the proton and the deuteron [24 
which implies that the ratio F^/F^ is not affected by higher twist contributions. We 



did not correct the data for target mass effects [25] so that these effectively are included 



in H(x). No higher twist, target mass or slow rescaling corrections |25|, |26| were applied 
to xFf e . 

The normalizations of the ZEUS, HI and NMC data sets were kept fixed to unity whereas 
those of E665, BCDMS, SLAC and CCFR were allowed to float within the quoted nor- 
malization uncertainties (7 parameters). There are thus in total 28 free parameters in the 
fit. 

The following cuts were made on the data: 

1. Discard data with x < s min = 0.001. This cut was introduced because the mass- 
less quark scheme used in this analysis is not suitable to describe low x structure 
functions where the contribution from charm is large. For the xF 3 data the cut was 
raised to x m i n = 0.1 to reduce the sensitivity to nuclear corrections at low x; 

2. Discard data with Q 2 < 3 GeV 2 . This cut reduces the sensitivity to QCD corrections 
beyond NLO. Because the QCD evolution of the ratio F^/F^ is small, we apply a 
lower Q 2 cut of 1 GeV 2 to these data; 

3. Discard data with W 2 = m 2 + Q 2 (l —x)/x < 7 GeV 2 (m p is the mass of the proton) 
to reduce the sensitivity to higher twist contributions and target mass effects which 
become important at high x and low Q 2 . 

The NLO QCD evolutions were calculated with the program Qcdnum P7fl; the \ 2 mini- 



mization and the calculation of the covariance matrices were done with Minuit 28 . 



3 Error propagation 

In this section we describe how the errors are defined in this analysis. The calculation 
of the statistical and systematic errors is based on linear error propagation. This implies 
that asymmetries in the input errors, if any, are ignored and that the output statistical 
and systematic errors are, per definition, symmetric. 
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In this analysis we have incorporated the effects of the point to point correlated experi- 
mental systematic errors in the model prediction for the structure functions. This model 
prediction, calculated at the kinematic point (xi,Qf), is defined as 

F l (p,s) = F? CD (p)(l + J2sAT) (5) 

3 

where F® CB (p) is the QCD prediction and A^J st is the relative systematic error on data 
point i stemming from source j.Q In Eq. (||) p denotes the set of parameters describing 
the input parton densities and s is a set of systematic parameters. Notice that the exper- 
iments, by giving central values and one standard deviation systematic errors, essentially 
provide 'measurements' of these systematic parameters: Sj — ± 1. 

Assuming that the Sj are uncorrelated and Gaussian distributed with zero mean and unit 
variance, the \ 2 can be written as 

where fi is the measured structure function and Afc is the statistical error with the point 
to point uncorrelated systematic errors added in quadrature. 

To propagate the statistical and systematic errors two Hessian matrices M and C were 
evaluated (with Minuit) at the minimum \ 2 : 



1 ay c _ i ay 

2 dpidpj ' 13 2 dpidsj 



The statistical covariance matrix of the fitted parameters is, as usual, given by 

ystat _ M -i 



A systematic covariance matrix of the fitted parameters can be defined as |29| 

V syst = M -l CC r M -l (g) 

The error on any function F of the parameters p (parton densities, structure functions, 
cross sections etc.) is then calculated using the standard formula for error propagation 

where V is the statistical, the systematic, or if the total error is to be calculated, the sum 
of both covariance matrices. 

The sources which contribute to the error on the results of the fit are: 



1. The statistical and systematic errors on the data. These were propagated as de- 
scribed above. The following systematic errors were taken into account in this anal- 
ysis: for the ZEUS data we included the 6 contributions as parameterized in 0. 

: If the systematic error is asymmetric the average of the upper and lower error is taken. 



5 



The SLAC and BCDMS systematic errors were taken to be those used in the QCD 
analysis of |24| . For the other data sets we included the systematic uncertainties as 
published (for HI and E665 only the normalization error and the total systematic 
error are available; the CCFR systematic errors are given in |30|| ). Apart from the 
7 normalization constants which were left free parameters in the fit, 50 independent 
sources of systematic error were propagated taking into account the correlations be- 
tween those of the NMC data sets. We assume, however, that the errors on F^j F^ 
are independent from those on F^ and F% . This has little effect on the total error 
estimate since the contribution from F 2 d / F% is small; 

2. The error on the strong coupling constant Aa s (M§) = 0.005, on the strange quark 
content of the proton AK S = 0.03 and on the nuclear corrections to F$, F%jF% 
and xF$ (AKd = 1, AKp e = 0.5). Also included is an assumed error on the charm 
threshold AQ 2 = 1 GeV 2 . To determine these uncertainties, the fits were repeated 
with the parameter lowered or raised by the error and the corresponding error band 
is defined as the envelope of the results from the central fit and these two additional 
fits; 

3. An additional error band ('analysis error') is defined as the envelope of the results 
from the central fit and several alternative fits described in Section |5|; 

4. The renormalization scale uncertainty. Two fits were performed with the renormal- 
ization scale set to fi 2 R = Q 2 /2 and fi 2 R = 2Q 2 while keeping the factorization scale 
fi 2 F fixed to Q 2 . The error band is defined as the envelope of the results from these 
fits and the central fit; 

5. The factorization scale uncertainty. Here the factorization scale was varied in the 
range Q 2 /2 < p? F < 2Q 2 while keeping the renormalization scale fi 2 R fixed to fi 2 F . 

Unless otherwise stated the errors shown in this paper are the quadratic sum of the 
contributions (l)-(3) listed above, that is, the scale errors are not included. The error 
contributions to the parton densities and several structure functions from each of the 
sources are given in Section |5] where also the scale uncertainties are shown. 

4 Results 

The fit yielded a good description of the data. Adding the statistical and systematic 
errors in quadrature a \ 2 = 1537 is obtained for 1578 data points and 28 free parameters. 
The x 2 values for each data set separately are listed in Table [TJ. The values of the fitted 
parameters are given in Table Table [l] shows that the \ 2 = 235/147 is rather high for 
the ZEUS F$ data. This is caused by 7 data points, more or less randomly distributed 
in the kinematic plane, which carry 84 units of x 2 - We have verified that removing these 
data points does not significantly change the results of the analysis. 

The good agreement between the F% data and the QCD fit can be seen from Fig. |2] where 
we show the x dependence of at fixed values of Q 2 . In Fig. |3] are plotted Ff and F 2 d from 
NMC, SLAC and BCDMS versus Q 2 for x > 0.1. Again, the agreement between the data 
and the fit is satisfactory. The full (dotted) curves in this plot show the QCD prediction 
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including (excluding) higher twist effects. The higher twist correction H(x) defined in 
Eq. (|D is plotted in Fig. [| (full curve). This correction is negative for x < 0.5 and becomes 



large and positive at high x. The same trend is observed in the analysis of MRST |3T 
(asterisks). We refer to Section [| for further investigations of higher twist effects. In the 
inset of Fig. |] we compare H(x) from this analysis with the result from [p^| , averaged 
over proton and deuteron. We observe that the higher twist contributions obtained in this 
analysis are substantially larger although part of the difference may be due to the fact 
that the authors of |24| have corrected the data for kinematical higher twist contributions 



(target mass effects) before extracting H(x). 

The Q 2 dependence of xF^ Fe for x > 0.1 is shown in Fig. |5|. There is good agreement 
between the data and the prediction from the QCD fit, corrected for nuclear effects as 
defined in Section |2|. Also the ratio F^j F% is well described as can be seen from Fig. [] 
where we plot the x dependence of F^/F^ = IF^/F^ — 1, averaged over Q 2 f!6 |. This 



figure shows that the error on the QCD prediction of F^/F^ rapidly increases above 
x = 0.4 indicating that the ratio of d to u quark densities, which is directly related to 
F2/F2, is not well constrained by the data at large x. 

We show in Fig. |7| the density x(d — u) as a function of x at Q 2 = 7.4 GeV 2 . There is 
reasonable agreement between the E866 data and the fit. For comparison we also plot 
the prediction from the CTEQ4 (CTEQ5 |l]]) parton distribution set which was obtained 
before (after) the E866 data became available. In Fig. |8] the strange quark density from 
the QCD fit is compared to the data from CCFR |23| (not included in the fit). The good 
agreement supports the description of x(s+s) with one fixed parameter K s (see Section g). 
The shaded bands in Fig. [8] show that the error on the strange quark density, which is 
dominated by that on K s , is adequately taken into account by assuming AK S = 0.03. 

In Fig. H are plotted the parton densities from the QCD fit versus x at Q 2 = 10 GeV 2 . 
Integrating the parton densities over x we find the following results for the momentum 
fractions carried by gluons and quarks at Q 2 = 4 GeV 2 : I g = f Q 001 xg dx = 0.394 ±0.018, 
I q = 0.594 ± 0.018 and I g + I q = 0.988 ± 0.002 (statistical and systematic errors added in 
quadrature) . The error on the last integral is much smaller than that on the first two since 
the total momentum fraction was constrained to be unity in the QCD fit. This clearly 
illustrates the importance of taking into account the correlations between the errors on 
the parton densities. Notice that 1 — I g — I q — 0.012 ± 0.002 is the small contribution to 
the momentum sum, at Q 2 = 4 GeV 2 , from the region x < 10~ 3 . 

Also shown in Fig. are the parton densities from CTEQ4. There is, within errors, good 
agreement. Notice, however, that xd v obtained in this analysis is slightly shifted towards 
larger values of x.f] Bodek and Yang || also obtain a harder xd density by modifying 
CTEQ4 (where the ratio d/u — > as x — > 1) such that 

d/u -> d'/u = d/u + Bx(l +x). (11) 

They find B = 0.10 ± 0.01 which implies that d/u — > 0.2 as x — > 1. 

The d/u ratio is shown in Fig. [1(| It is seen that the result from the QCD fit with B = 
(full curve) for x < 0.75 is close to the modified CTEQ4 distribution with B = 0.1 (dashed 
curve) . If we parameterize the d quark density according to Eq. (|ll]) and leave B a free 
parameter in the fit we obtain B = —0.02 ± 0.01 (statistical error), close to zero. In any 

2 A similar shift is seen in the most recent parton distribution sets CTEQ5 and MRST 0]. 
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case, the large error band shown in Fig. 10 clearly indicates that the exact behavior of 
d/u at large x is presently not well constrained by the data. This ratio might go to zero 
(CTEQ), to a constant (Bodek and Yang) or may even diverge (this analysis) as x —> 1. 



5 Systematic checks 

To check the stability of the results we have performed a series of alternative fits: 



1. Use the quadratic sum of the statistical and systematic errors in the definition of 
the x 2 instead of taking statistical errors only; 

2. Release the momentum sum rule constraint. In this fit the total momentum fraction 
carried by quarks and gluons is found to be 1.03, close to unity. We did not perform a 
full error analysis for this fit but the uncertainty on the momentum sum is probably 
about the same as that given in Section |] for the gluon and quark momentum 
fractions (±0.02); 

3. Set the value of the input scale to Qq = 7 instead of to 4 GeV 2 ; 

4. Lower the Q 2 cut from 3 to 2 GeV 2 and the W 2 cut from 7 to 5 GeV 2 ; 

5. Raise the Q 2 cut from 3 to 4 GeV 2 and the W 2 cut from 7 to 10 GeV 2 ; 

6. Remove the lower cut of x m \ n = 0.1 on the xF 3 data. 

7. Fix the normalization of the CCFR xF 3 data (to iV = 1.008, see Table |l|) and let 
the nuclear correction on xF£ Fe float by leaving Kp e a free parameter in the fit. 
This results in Kp e = 0.8, consistent with the input value of 1.0 ± 0.5; 

8. Use the parameterization of Seligman |30] (see Section [!]) instead of that from Es- 
kola et al. to correct the xF^ Fe data for nuclear effects; 

9. Leave the normalization of the NMC data free (one parameter for the 8 NMC data 
sets) so that all fixed target data are re-normalized with respect to HERA. This 
results in A^nmc — 1-007 while the normalizations of the other data sets changed by 
less than 0.01, compared to the central fit; 

10. Parameterize the d v density as given in Eq. ( |TT1) and leave the parameter B free in 
the fit. As already mentioned in Section |] we find B = —0.02 ± 0.01 (statistical 
error) . 



All the fits listed in items (1)-(10) above yielded good descriptions of the data with a x 2 
per degree of freedom close to that of the central fit. To quantify the stability of the QCD 
analysis we define an 'analysis error band' as the envelope of the results from the central 
fit and these 10 alternative fits. 

In Table ^ the contributions to the relative error on several parton densities and structure 
functions are given for each source listed in Section £| These error contributions are 
calculated at Q 2 = 10 GeV 2 and are averaged over x. In Fig. Ill] we show as a function 
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of x at Q 2 = 10 GeV 2 the contributions to the relative error Ag/g (left-hand plots) and 
AE/E (right-hand plots). It is seen that the experimental errors shown in Fig. |lTJa are 
roughly equal to the input errors which, in Fig. |TT|b, are calculated as the quadratic sum 
of the contributions listed in item (2) of Section |3|. Compared to these two contributions 
the analysis error band is small (Fig. |TT|c). This illustrates the stability of the fit with 
respect to the ten alternatives listed above. The factorization and renormalization scale 
uncertainties added in quadrature are shown in Fig. |TT|d. For the gluon distribution the 
scale error is roughly equal to the other error contributions. For the quarks however 
the scale error, which is dominated by the factorization scale uncertainty, is the largest 
contribution to the total error. 

To check the higher twist contributions obtained from the central fit, we have determined 
these in an alternative way as follows. First, we performed a fit with the Q 2 and W 2 cuts 
raised to 10 and 35 GeV 2 respectively. This restricts the kinematic region to a domain 
where higher twist effects are small (< 4% as estimated from the results of the central fit) 
so that the Q 2 evolution of the data can be described by perturbative QCD only. Keeping 
H(x) to zero in the fit, a good description of the data was obtained with a x 2 = 824 for 
901 data points and 21 free parameters.^ The parton densities obtained from this fit are, 
within errors, consistent with the central result. 

Second, we used the parton densities obtained above to calculate the QCD predictions 
for F% and F 2 d down to the standard cuts of Q 2 = 3 and W 2 = 7 GeV 2 . Keeping F 2 QCD 
in Eq. (Q) fixed to these predictions we determined H(x) from a fit to the data on F\ 
and F$. The result is compatible with that obtained from the central fit as shown by the 
dashed curves in Fig. f|. We conclude from this investigation that no significant bias is 
introduced by including, in the central fit, higher twist effects as parameterized in Eq. (^). 

Finally, we performed a fit (with the standard cuts on the data) where the momentum 
sum rule constraint was removed and where in addition the higher twist contribution 
H{x) was kept fixed to zero. This resulted in a considerably worse x 2 = 1918 for 1578 
data points and 24 free parameters. Furthermore the total momentum fraction carried 
by quarks and gluons was found to be 1.20. This is so because a large gluon density 
is favored by the large slopes of the F 2 data at intermediate x and low Q 2 (see Fig. [|). 
We reject this fit as an alternative to the central result because of the large x 2 an d the 
violation of the momentum sum rule. 



6 Comparisons 

In this section we calculate a few predictions based on the parton densities obtained from 
the QCD analysis and compare these to deep inelastic scattering data which are not 
included in the fit (except xF^ Fe at large x). 

First, we investigate nuclear effects on the u-Fe structure functions by calculating the ratio 
of the CCFR data to the predictions from the QCD fit of neutrino scattering on a free 
nucleon. We find that the ratios xF^ Fe / xF^ N and F^ Fe jF^ N do not, within errors, depend 
on Q 2 in accordance with measurements of nuclear effects in charged lepton scattering 

3 In addition to the five higher twist parameters the normalizations of the SLAC data sets were kept 
fixed because these data were almost entirely removed by the larger cuts. 
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where also no significant Q 2 dependence has been observed [32]. We have therefore fitted 
the xF 3 data to the parameterization 



xF% (%i,Qj) — Ri xF^ {xiiQj) 

where xF^ N was kept fixed to the QCD prediction and Ri = xF^ Fe / xF^ N \xi) was left 
a free parameter for each bin of x. In the fit only statistical errors on the data were 
taken into account. The upper (lower) systematic error on Ri was estimated by offsetting 
the data with each source of systematic error in turn, repeating the fit and adding the 
deviations above (below) the central value in quadrature. Likewise the ratio F% Fe : /K 



2 



was calculated for each bin of x. In the fits only data were included above the charm 
threshold Q 2 C = 4 GeV 2 . 

In Fig. [12| the results are plotted as functions of x for xFg /xFg (top) and F 2 I/Fe /F 2 I/Ar 
(bottom). Both these ratios exhibit the typical x dependence of nuclear effects including a 
hint for the rise at large x due to Fermi motion in the nucleus. It is clear that nuclear effects 
in neutrino scattering on heavy targets are present and that corrections must be applied 



if these data are included in QCD fits. Also shown in Fig. [12| are two parameterizations 
of nuclear corrections from Eskola et al. (full curves) 4 and from Seligman |$0| (dashed 
curves). Both these parameterizations are obtained from fits to data on nuclear effects in 
charged lepton scattering. To show the contribution to the errors on xF^ Fe / xF^ N and 
F 2 I/Fe /F 2 l/Ar from the QCD fit we have chosen to draw these as the shaded bands around 



the full curves. We conclude from Fig. 12 that, in the kinematic range explored and 



within the present errors, there is no strong indication for nuclear effects being different 
in neutrino or charged lepton deep inelastic scattering. 

Next, we investigate to which extent the massless quark scheme used in this analysis is 
able to describe the charm contribution F% to the inclusive F 2 structure function. In 
Fig. [13] we show the recent data from ZEUS on F 2 |3|| versus Q 2 for x > 5 x 10~ 4 . The 



full curves with the shaded error bands correspond to the prediction of _F 2 from the QCD 
fit. Because these predictions are calculated with the assumption that the charmed quark 
is massless, the curves extrapolate to _F 2 C = at the charm threshold Q 2 = 4 GeV 2 in 
disagreement with the data. For comparison we also plot in Fig. [T3| the result from the 



ZEUS QCD analysis of [34| where charm mass effects were taken into account using the 
heavy quark three flavor scheme as defined in JlO| (dashed curves). This heavy quark 
scheme gives a good description of F 2 C over the full Q 2 range down to the lowest measured 
value of x = 5 x 10^ 5 , see [fffi| . However, it is seen from Fig. [13| that above Q 2 ~ 10 GeV 2 
also the massless quark scheme describes the data well. This supports the assumption 
made in the remainder of this section that, at least for x > 10~ 3 , quark mass effects do 
not spoil the QCD extrapolations to large Q 2 . 

To compare with the recent e + p data from ZEUS M, [5| we have used the parton densities 
from this analysis to calculate the Born cross sections for the processes e + p — > e + X (NC) 
and e + p — > vX (CC). For unpolarized positrons and protons these cross sections are, at 
HERA energies, to good approximation given by 

d 2 a NC (e + p) _ 27ra 2 r NC 2 NC NC -, 

dxdQ 2 ~ xQ* [Y+t2 V Y - Xt * J {12) 

vN 



4 Plotted are the nuclear corrections for iron on xu v + xd v oc xFg N (see Section |[) and xT, oc F%* 



10 



d 2 a cc (e+p) G 2 F M* 

dxdQ 2 4^ (Q 2 + M' w 



K \v, F cc - v 2 F cc - Y tF cc ^ 



Here a is the fine structure constant, Gp the Fermi constant and M^y the mass of the W 
boson. In Eq. (|12"D , y = Q 2 /xs and Y± — 1 ± (1 — y) 2 with s the square of the e + p center 
of mass energy. The structure functions F2, Fl and xF% were calculated as functions of 
x and Q 2 in NLO using the expression given in [35| for Fl. In the NC case contributions 
from Z° exchange and 7-Z interference were taken into account. We refer to Q and || 
and references therein for further details concerning the calculation of the NC and CC 
structure functions and for the input values of the electroweak parameters. We mention 
here that in leading order QCD the CC cross section can be written as |36| 



d 2 a cc {e + p) _ G 2 F Myy 

f2 \2 



W 



dxdQ 2 ~~ 2ttx (Q 2 + M, 



w, 



[(1 — y) 2 [xd + xs) + xu + xc\ (13) 



from which it is seen that this cross section is predominantly determined by the down 
quark density. 

In Fig. [14| we plot the ZEUS data on the NC e + p single differential cross sections da/dQ 2 , 
da/dx and da/dy for Q 2 > 400 GeV 2 together with the predictions from the QCD fit 
(full curves). Here both the data and the predictions are divided by the cross sections 
calculated from the CTEQ4D parton distribution set. This figure shows that the present 
analysis and CTEQ4 give a good description of the measurements although there is a 
slight tendency (~ 5%) of the QCD predictions to be below the data at low Q 2 , low x 
and large y. These are the regions where the calculation of the cross sections might be 
sensitive to higher order QCD corrections on Fl |37] • 



The ZEUS results on the CC e + p single differential cross sections for Q 2 > 200 GeV 2 
are shown in Fig. |15| (again divided by the CTEQ4 predictions). The agreement between 
the data and the predictions from the QCD fit (full curves) is excellent. The larger cross 
section at large x and Q 2 , compared to CTEQ4, is due to the harder d v distribution 
obtained in this analysis, see Fig. [t| We have verified that using the more recent parton 
distribution sets CTEQ5 and MRST, also with a harder d density, yields results close to 
those presented here. 

We conclude from Figs. Ill and ITS] that, within the present accuracy, no significant devi- 



ations can be observed between the data and the Standard Model predictions for Q 2 up 
to about 3 x 10 4 GeV 2 . 



7 Summary 

The parton momentum density distributions in the proton were determined from a NLO 
QCD analysis of structure functions measured at HERA and by fixed target experiments, 
together with measurements of the difference x(d — u). The data included in the fit cover 
a kinematic range of 0.001 < x < 0.75 , 3 < Q 2 < 5000 GeV 2 and W 2 > 7 GeV 2 . The 
structure functions measured in nuclear targets were corrected for nuclear effects. Higher 
twist contributions to the proton and deuteron F 2 structure functions were taken into 
account phenomenologically. The fit was performed in the variable flavor number scheme 
where the charm and bottom quarks are assumed to be massless. It turns out that this 
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scheme gives for x > 10 3 and Q 2 > 10 GeV 2 a good description of the charm structure 
functions measured at HERA. 

The uncertainties in the parton densities, structure functions and related cross sections 
were estimated from the experimental statistical errors, from the 57 independent sources 
of systematic uncertainty and from the errors on the input parameters of the fit. An 
additional analysis error is defined as the envelope of the results from the central fit and 
10 alternative fits. This error gives a small contribution to the total error. Also estimated 
are the renormalization and the factorization scale uncertainties. 

The fit yields a good description of the data. The higher twist contributions to F% and 
F$ are found to be significant over the whole x range covered in this analysis. The ratio 
of down to up quark densities at large x > 0.4 is ill constrained by the data. As a 
consequence we find little sensitivity to a recently proposed parameterization of the d 
quark density aimed at modifying d/u at high x. 

Nuclear effects in neutrino scattering on iron were investigated by comparing measure- 
ments of xF^ Fe and with the QCD predictions of neutrino scattering on a free nu- 
cleon. We do not, within the present errors, find a significant difference between nuclear 
effects in neutrino and charged lepton scattering. 

From the evolved parton densities the e + p NC and CC structure functions and Born 
cross sections were calculated. There is good agreement between these Standard Model 
predictions and the measurements by ZEUS of the e + p NC and CC cross sections at large 
x and Q 2 . 

The results of this analysis are stored on a computer readable file which is available from 



the web site prttp : //www.nikhef .nl/user/h24/qcdnuiii| . This file contains the statistical 
and systematic covariance matrices and, as functions of x and Q 2 , the parton densities 
as well as the derivatives of these densities to the fitted parameters. Also stored are the 
results of the systematic fits given in Section § item (2)-(5) as well as the analysis error 
band defined in Section [5]. A Fortran program |38| is provided to read the file and it 
contains tools to propagate the errors to any function of the parton densities. 
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Dataset 


F 


x 2 


points 


X 2 /point 


norm. 


SLAC P 




61.3 


56 


1.10 


0.988 


SLAC D 


pi 


52.3 


57 


0.90 


0.984 


BCDMS P 




155.2 


177 


0.88 


0.979 


BCDMS D 


Fi 


159.3 


159 


1.00 


1.000 


NMC E090 P 


n 


46.0 


44 


1.04 




NMC E120 P 




63.2 


53 


1.19 




NMC E200 P 




75.4 


64 


1.18 




NMC E280 P 




65.2 


72 


0.91 




NMC E090 D 


Fi 


47.2 


44 


1.07 




NMC E120 D 


Fi 


43.5 


53 


0.82 




NMC E200 D 


Fi 


47.6 


64 


0.74 




NMC E280 D 


Fi 


51.4 


72 


0.71 




E665 P 


Fi 


53.4 


41 


1.30 


1.017 


E665 D 


Fi 


44.6 


41 


1.09 


1.001 


ZEUS 


Fi 


235.1 


147 


1.60 




HI 


Fi 


98.2 


150 


0.65 




NMC 


Fi/Fi 


189.1 


205 


0.92 




CCFR 


xF vFe 


33.2 


68 


0.49 


1.008 


E866 


x(d — u) 


16.9 


11 


1.54 




Total 




1537.2 


1578 


0.97 





Table 1: The contribution to the total x 2 from each data set included in the QCD fit. The 
X 2 values are calculated with the statistical and systematic errors added in quadrature. 
The last column gives the values of the normalization constants left free in the fit. 



Parameter xg xS xA xu l 



A 
5 
V 

7 



1.32 
-0.26 ±0.03 

5.19 ±1.53 
-0.52 ±1.42 



0.81 ±0.03 
-0.15±0.01 

3.96 ±0.19 
-1.32±0.06 



0.31 ±0.10 
0.57±0.09 
7.47 ±1.00 



2.72 

0.62 ±0.02 
3.89 ±0.02 
2.79 ±0.31 



C^/^j 

1.98 

0.65 ±0.02 
3.07±0.18 
-0.82 ±0.12 



h hi fi2 hs h 4 

-0.12±0.06 —1.9 ±1.1 12.7±5.7 -45.0 ±12.7 55.2±9.9 



Table 2: The values of the parameters obtained from the QCD fit. The parameters 
are defined in Section |2|. The errors given are statistical only. The full statistical and 
systematic error matrices can be accessed as described in Section [7]. 
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/ 


stat. 


syst. 


a s 


K a 


K d 


K Fe 


Ql 


anal. 


total 


A 




xg 


3.1 


7.1 


12.1 


0.1 


1.3 


1.8 


3.3 


4.7 


16.1 


10.1 


8.5 


xH 


0.7 


3.9 


1.2 


1.2 


0.3 


0.8 


2.5 


1.6 


5.5 


0.6 


5.8 


xS 


1.1 


5.4 


1.4 


1.6 


0.6 


1.7 


3.5 


2.5 


8.1 


0.9 


7.9 




0.9 


3.4 


1.5 


0.0 


1.3 


3.6 


0.1 


1.0 


5.7 


2.1 


1.4 




1.7 


5.8 


1.9 


0.0 


2.3 


6.4 


0.3 


4.2 


11.0 


2.6 


1.4 


Fl 


0.3 


3.3 


0.1 


0.0 


0.1 


0.2 


0.1 


0.8 


3.4 


0.2 


0.3 


Fi 


0.3 


3.3 


0.1 


0.0 


0.4 


0.3 


0.2 


0.8 


3.5 


0.2 


0.3 


Ft m 


0.1 


0.4 


0.0 


0.0 


0.3 


0.1 


0.0 


0.1 


0.6 


0.1 


0.0 


xF vFe 


0.9 


3.6 


1.5 


0.0 


1.1 


3.6 


0.1 


3.5 


7.2 


2.4 


1.0 


xA 


9.7 


21.7 


1.9 


0.2 


29.0 


9.5 


5.2 


14.4 


45.5 


10.9 


4.1 



Table 3: The relative errors Af/f (in percent) on the parton densities and several 
structure functions, averaged over x. The errors are calculated at Q 2 = 10 GeV 2 . The 
columns give the individual contributions from the sources described in Section |3|. The 
column marked 'total' gives the total error excluding the scale errors. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 X 



Figure 1: The ratio of nuclear structure functions to those of a free nucleon versus x: 
(a) F$/Fg from @; (b) xFg Fe /xF£ N from |0|. The shaded bands correspond to the 
uncertainties on the ratios assumed in this analysis. Also indicated are the x ranges of 
the F 2 and xF 3 data included in the QCD fit. 
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Figure 2: The proton F 2 structure function versus x at fixed values of Q 2 from ZEUS 
HI 0, E665 @, NMC @, BCDMS @ and SLAC pf. Only shown are the c 
included in the QCD analysis. The full curves show the result from the QCD fit. 



clarity the constants given in brackets are added to F 2 . The ZEUS, HI and E665 data 
points are plotted with open and solid symbols for alternating Q 2 bins. 
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Figure 3: The structure functions F% (left) and (right) from fixed target experi- 
ments [TJ, |T3| versus Q 2 for x > 0.1, Q 2 > 3 and W 2 > 7 GeV 2 . The full (dotted) 



curves correspond to the QCD fit results including (excluding) higher twist contributions. 
The F$ curves include the correction for nuclear effects in deuterium described in the 
text. The x values and multiplication factors (in brackets) given in the left-hand plot 
apply to both F$ and F$. 
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Figure 4: The higher twist correction H(x) as defined in Section |2|. The full curve corre- 
sponds to the result of the QCD analysis. The quadratic sum of the error contributions 
(1) and (2) given Section ^| is drawn as the shaded band around the curve. The dashed 
curve shows H{x) obtained from an alternative fit described in Section [|. Also plotted 
are the results from an analysis of MRST |31| (asterisks). In the inset H(x) from this 
analysis is compared to the result from the QCD fit of |23| to the SLAC and BCDMS 
data. 
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Figure 5: The xF^ Fe structure function from CCFR Jl7| versus Q 2 at fixed values of 
x > 0.1. The full curve shows the QCD prediction corrected for nuclear effects on iron, 
described in the text. The shaded bands indicate the error on the QCD fit. For clarity 
xF 3 is multiplied by the factors indicated in brackets. 



21 



v 2 /v 2 

1.1 



0.9 h 

0.8 

0.7 h 
0.6 

0.5 - 

0.4 - 



0.05 


-0.05 h 



NMC 
QCD fit 



_i i i i_ 



_L 



Data - Fit 



_i i i ' ' ' 




10 



10 



X 



Figure 6: Top: The ratio F^/Ff versus x from NMC, averaged over Q 2 JT^]. The full 
curve shows the prediction from the QCD fit, corrected for nuclear effects in deuterium 
as described in the text. Bottom: The difference between the data and the QCD fit. The 
shaded band shows the error on the QCD prediction. 
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Figure 7: The difference x(d — u) from E866 versus x at Q 2 = 7.4 GeV 2 . The 
full curve shows the result from the QCD fit with the error drawn as the shaded band 
around the curve. The dashed (dashed-dotted) curve is the prediction from CTEQ4 || 
(CTEQ5 0). 
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x(s + s) 




Figure 8: The strange quark density x(s + s) versus x at three values of Q 2 . The full 
circles show the data from CCFR [^3J. The curves correspond to the result from the QCD 
fit. The error on the QCD prediction is indicated by the shaded bands. 
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Figure 9: The parton momentum densities xg, xS (both divided by a factor of 20), 
xu v and xd v versus x at Q 2 = 10 GeV 2 . The full curves show the results from the QCD 
fit with the errors drawn as shaded bands. The dashed curves are from the CTEQ4 || 
parton distribution set. 
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Figure 10: The ratio (d + d)/(u + u) versus x at Q 2 = 10 GeV 2 from the QCD fit (full 
curve) and CTEQ4 || (dashed-dotted curve). The shaded band shows the error on the 
QCD fit. The dashed curve corresponds to the CTEQ4 prediction with a modified down 
quark density as described in the text. Indicated is the x range covered by the data 
included in the QCD fit. 
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Figure 11: The four contributions to the relative error on the gluon density Ag/g (left 
hand plots) and the quark density AS/S (right hand plots) versus x at Q 2 = 10 GeV 2 . 
The error contributions are defined in the text. 
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Figure 12: The ratios xF» Fe /xF» N (top) and F% Fe /F» N (bottom) of the CCFR u-Fe 
data |T7| to the prediction from the QCD fit for neutrino scattering on a free nucleon, ob- 
tained from a fit described in the text (solid dots). The errors correspond to the statistical 
and systematic errors on the data added in quadrature. Also shown are parameteriza- 
tions of nuclear effects in charged lepton scattering on iron from Eskola et al. (full 
curves) and Seligrnan |JU| (dashed curves). The error contributions from the QCD fit to 
xF^ Fe /xF% N and F% Fe / F% N are drawn as the shaded bands around the full curves. The 
x range of the xF 3 data included in the QCD analysis is indicated in the upper plot. 
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Figure 13: The charm structure function F 2 C measured by ZEUS [33] versus Q 2 at fixed 
values of x > 5 x 10~ 4 . The full curves with shaded error bands show the predictions from 
the QCD fit. These curves extrapolate to F| = at the charm threshold Q 2 C = 4 GeV 2 . 
The dashed curves correspond to the predictions from the ZEUS QCD fit |34| where charm 
mass effects were taken into account. For clarity F 2 C is multiplied by the factors indicated 
in brackets. 
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Figure 14: The e + p NC single differential cross sections da/dQ 2 (top), da/dx (middle) 
and da/dy (bottom) measured by ZEUS ||, divided by the QCD prediction calculated 
with the CTEQ4D || parton distribution set. The full curves with shaded error bands 
correspond to the results of the QCD fit. The inset in the top plot shows da/dQ 2 for 
Q 2 <2x 10 4 GeV 2 . 
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Figure 15: The e + p CC single differential cross sections da/dQ 2 (top), da/dx (middle) 
and da/dy (bottom) measured by ZEUS f|, divided by the QCD prediction calculated 
with the CTEQ4D || parton distribution set. The full curves with shaded error bands 
correspond to the results of the QCD fit. 
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